Search for: All records

Monitoring the activity of subglacial volcanoes along the Aleutian Arc in Alaska is important to the safety of local populations, as well as air traffic flying through the region. However, observations of volcanic unrest are limited by accessibility and resources, particularly at glacier-covered systems, making investigations of their stability challenging. Westdahl Peak, a subglacial volcano on Unimak Island in the Aleutian Arc has experienced significant unrest and uplift since its most recent VEI three eruption in 1991-1992. Given the magnitude of observed uplift, previous investigations suggested the potential for eruption by 2010, but no such event has occurred. One hypothesis to explain this prolonged unrest is that the 1-km thick glacier may increase the stability of the magma system. However, the impact of ice caps and glaciers on the short-term stability of volcanoes is not well understood. In this study, thermomechanical finite element models are used to evaluate how the stability of a glaciated volcano is impacted by variations in ice cap thickness, magma chamber depth, geometry, magma flux rate, and seasonal changes in ice cover thickness. Our numerical experiments indicate that the presence of an ice cap (1–3 km thick) increases the average repose interval for a magma system. Among models with different magma chamber geometries, depths, and flux rates, the greatest increases in repose interval are observed in prolate systems where the increase is up to 57% for a chamber located at 5 km-depth. Spherical and oblate also experience smaller, yet significant, increases in repose interval. Additionally, the percentage increase in repose interval is not impacted by variations in magma flux rate for a given ice cap thickness and magma chamber geometry. However, flux rates do influence the timing of eruptions when the system is experiencing seasonal variations in ice thickness. Our results show that systems with low flux rates are more likely to fail when the ice thickness is at its lowest. The numerical estimates further suggest that the ice cap on Westdahl Peak, which is ∼1 km, may slightly increase the stability of the magma system. In general, given flux rates and magma chamber geometries estimated for the Westdahl system, the repose interval can increase by ∼7 years due to the Westdahl glacier. This increase is small on a geologic scale but is significant on human time scales and the impact of glaciers must be considered in future forecasting efforts. 

Ensemble based data assimilation approaches, such as the Ensemble Kalman Filter (EnKF), have been widely and successfully implemented to combine observations with dynamic forecast models. In this study the EnKF is adapted to assimilate ground deformation observations from interferometric synthetic-aperture radar (InSAR) and GPS into thermomechanical finite element models (FEM) to evaluate volcanic unrest. Two eruption hindcasts are investigated: the 2008 eruption of Okmok volcano, Alaska and the 2018 eruption of Sierra Negra volcano, Galápagos, Ecuador. At Okmok, EnKF forecasts tensile failure and the lateral movement of the magma from a central pressure source in the lead up to its 2008 eruption indicating potential for diking. Alternatively, at Sierra Negra, the EnKF forecasts significant shear failure coincident with a Mw 5.4 earthquake that preceded the 2018 eruption. These successful hindcasts highlight the flexibility and potential of the volcano EnKF approach for near real time monitoring and hazard assessment at active volcanoes worldwide. 

Understanding the stress evolution of extinct volcanoes can improve efforts to forecast flank eruptions on active systems. Field, petrographic, and seismic data are combined with numerical modeling to investigate the paleo‐stress field of New Zealand's Akaroa Volcano, or Akaroa Volcanic Complex. Field mapping identifies 86 radially oriented dikes and seven lava domes found only within a narrow elevation range along Akaroa's erosional crater rim. These observations suggest that crater rim dike emplacement resulted from lateral deflection of vertically ascending intrusions from a centralized magma source, which in turn may have facilitated formation of the lava domes, as well as two scoria cones. We postulate that dike deflection occurred along a stress barrier, as neither a compositional change nor structural boundary are present. We use a finite element model (FEM) simulating Akaroa to test how different factors may have influenced the system's stress state and dike geometry. Elastic, non‐flexural (“roller”) model configurations containing a large, oblate, and shallow magma chamber produce stress barriers most conducive to radial dike emplacement along Akaroa's crater rim. These configurations also simulate rapid edifice construction above a preexisting lithospheric “bulge.” Conversely, simulating flexural stresses exerted on the lithosphere by Akaroa's large mass hinder rather than promote radial dike emplacement. Temperature‐dependent viscoelastic relaxation promotes gradual increases in stress barrier elevation, though this effect is strongly dependent on magma chamber parameters. These results suggest that Akaroa was constructed rapidly (within ∼100 kyr) prior to crater rim dike emplacement, which occurred throughout the volcano's remaining active lifespan.